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The shape of hadronic distribution amplitudes (DAs) is a critical issue for the 
perturbative QCD of hard exclusive processes. Recent CLEO data on 77* — » 7T° 
form factor clearly favor a pion DA close to the asymptotic form. We argue that 
QCD sum rules for the moments of the pion DA ip-^(x) are unreliable, so that 
the humpy shape of tp^ix) obtained by Chernyak and Zhitnitsky is a result of 
model assumptions rather than an unambigous consequence of QCD sum rules. 
This conclusion is also supported by a direct QCD sum rule calculation of the 
77* — » 7T° form factor which gives a result very close to the CLEO data. 



1 Introduction 

In this talk, I discuss some general features of QCD sum rule applications to 
hadronic wave functions and form factors using as examples the pion distribu- 
tion amplitude ^^(x) and transition form factor for the process 7*7* — > 7r° in 
which two virtual photons produce a neutral pion. This process provides an 
exceptional opportunity to test QCD predictions for exclusive processes. In 
the lowest order of perturbative QCD, its asymptotic behaviour is due to the 
subprocess 7*(<?i) + 7* (92) — » q(xp) + q(%p) with x (x) being the fraction of 
the pion momentum p carried by the quark produced at the q\ (92) photon 
vertex. The relevant diagram resembles the handbag diagram of DIS with 
the pion distribution amplitude (DA) ^>tt1.x) instead of parton densities. The 
asymptotic PQCD prediction is given byld (q 2 = —q 2 , q\ = —Q 2 , x = 1 — x): 

F 7 ,^ [q,Q)- 3 ^ xQ2+ _^ dx ^ 3 ^ dx _ 5Q2 I. (l) 

Experimentally, thejnost important situation is when one of the photons is 
almost real q 2 ~ Era. In this case, necessary nonperturbative information 
is accumulated in the same integral I (see Eq.(|l|)) which appears^ixL.the one- 
gluon-exchange diagram for the pion electromagnetic form factor Eat 
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of / is sensitive to the shape of the pion DA ip v (x) . Using the 
form tf%(x) = 6Uxx gives / = 3 and F™ {Q 2 ) = AttU/Q 2 . 



The va 

asymptotic-. ........ r , „ 

If one takes the CZ form til ip^ z (x) — 30/ w xx(l — 2x) 2 , then 1 = 5, and this 
difference can be used for experimental discrimination betwaen-the two forms. 
One-loop radiative QCD corrections to Eq.(jl]) are known IStM and they are 
under control. Clearly, the asymptotic l/Q 2 -behaviour cannot be true in the 
low-Q 2 region, since the Q 2 = limit of F 1 yt„o(Q 2 ) is known to be finite 
and normalized by the ir° — > 77 decay rate. From the axial anomaly lij, 
F 77 » 7r o(0) = l/nfir- Brodsky and Lepage El proposed a simple interpolation 
7tUF^^ (Q 2 ) = 1/(1 + Q 2 /4tt 2 /2) between the Q 2 = value 1/tt/ w and the 
leading-twist PQCD behavior 47r/ w /Q 2 with normalization corresponding to 
the asymptotic DA. Note that the mass scale = 47r 2 / 2 w Q.67 GeV 2 in 
this monopole formula is close to m 2 . Recent experimental datat! from CLEO 
are below the BL-curve and are by almost a factor of 2 lower than the value 
for the CZ wave function. This result apparently excludes the CZ DA and 
suggests that the pion DA may be even narrower than ip^(x). Since the CZ 
model is often perceived as a direct consequence from QCD sum rules, the 
experimental evidence in favor of a narrow DA may be treated as a failure of 
the QCD sum rule approach. One should remember, however, that accuracy 
of QCD sum rules strongly depends on the specific hadronic characteristics to 
which the sum rule technique is applied. Long ago, in papersEil written with S. 
Mikhailov, we argued that CZ sum rules are very unreliable, with the results 
strongly depending on the assumptions about the size of higher terms in the 
operator product expansion (OPE). 

2 QCD sum rules 

QCD sum rulesEl are based on quark-hadron duality, i.e., possibility to describe 
the same object in terms of either quark or hadronic fields. To calculate f v , 
consider the p^p^-part of the correlator of two axial currents: 

iP» = i J e ^(o|T(j+(^)i 5 -(o))|o)rf^ = p^n 2 (p 2 )-^ I/ n 1 ( P 2 ). (2) 

The dispersion relation represents Ll2(p 2 ) as an integral over hadronic spectrum 
2 1 f 00 p hadTOn (s) 



n 2 (p ) = - / ^-ds + "subtractions" (3) 

K Jo S ~P 

with the spectral density p hadron (s) determined by projections of the axial 
current onto hadronic states ((01^(0)1^; P) = if^P^, etc.): 

p hadro„ (s) _ nftSia - ml) + Trftjis - m 2 Al ) + "higher states" (4) 
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j-yexp K MeV in our normalization). On the other hand, when the probing 
virtuality is negative and large, one can use the OPE 



n 2 ( P 2 ) = nrV) + 4kgg) + ^a s (q q f + ... (5) 



where n^'fjJ 2 ) = Ii2 Ualk (p 2 ) is the perturbative version of I±2(p 2 ) given by 
a sum of PQCD Feynman diagrams while the condensate terms (GG), (qq), 
etc., (with calculable coefficients A, B,etc. ) describe/parameterize the non- 
trivial structure of the QCD vacuum. The quark amplitude II quark (p 2 ), can 
also be written in the dispersion representation (0), with p(s) substituted by its 
perturbative analogue /o quark (s) = 4- (l + ^ + .. •) (quark masses neglected). 
Hence, the condensate terms describe the difference between the quark and 
hadron spectra. Treating the condensate values as known, one can try to 
construct a model for the hadronic spectrum. The simplest model is to ap- 
proximate all the higher resonances including the Ai by the quark spectral 
density starting at some effective threshold sq: 

Neglecting the pion mass and using the standard values for the condensates 
(GG), (qq) 2 , one should adjust sq to get an (almost) constant result for the 
rhs of the SVZ-borelized version of the sum rule 

,2 _ 1 f _quark M p -s/M* . , j^gg) , 176 Tra s {qq} 2 

U ~nJ 9 [) + + 81 Af 4 + -''' Uj 

The magnitude of /„■ extracted in this way, is very close to its experimental 
value /° xp « 130 M el/. Changing the values of the condensates, one would get 
the best M 2 -stability for a different so, and the resulting value of /„. would 
also change. Correlation between the fitted values of f v and sq is manifest in 
the M 2 — > oo limit of the sum rule 
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so 



# = -/ p^' k ( s )d s , (8) 



giving a local duality relation which states that two densities p quark (s) and 
p hadron (s) give the same result if one integrates them over the appropriate 
duality interval sq. The role of the condensates was to determine the size of the 
duality interval sq, but after it was fixed, one can write down the relation 
which does not involve the condensates. In the lowest order, Pg Uark (s) = 1/47T, 
which gives so = 47r 2 / 2 . Note, that this is exactly the combination which 
appeared in the Brodsky-Lepage interpolation formula. 
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3 CZ sum rules and pion DA 

Chernyak and A. Zhitnitskyi proposed to use QCD sum rules for calculating 
next moments (£ N ) (where £ = 2x — 1) of the pion DA. They extracted (£ 2 ) 
and (£ 4 ) from the relevant SR 

f 2 ( , N) 3A/ 2 (1-e-WM') g s (GG) 16 7ra.(g g ) a 

M * ' 4tt 2 (A + l)(A + 3) 12ttM 2 + 81 M 4 ^ ii + 4JV J W 

precisely in the same way as the f„ value. Note that the scale determining the 
magnitude of the hadronic parameters is settled by the ratios of the condensate 
contributions to the perturbative term. If the condensate contributions in the 
CZ sum rule (|^) would have the same A-behavior as the perturbative term, 
then the A-dependence of (£ N ) would be determined by the overall factor 
3/ (A + 1)(A + 3) and the resulting wave function <p% s (x) = 6f w x(l — x) would 
coincide with the asymptotic form. However, the ratios of the (qq) and (GG)- 
corrections to the perturbative term in Eq. (^) are growing functions of A. In 
particular, in the (qq) case, the above mentioned ratio for A = 2 is by factor 
95/11 larger than that in the A = case. For A = 4 the enhancement factor 
equals 315/11. As a result, the effective vacuum scales of (mass) 2 dimension 
are by factors (95/11) 1 / 3 « 2.1 and (315/11) 1 / 3 w 3.1 larger than that for the 
A = case. Approximately the same factors (5 1 / 2 « 2.2 and (35/3) 1 / 2 w 3.4) 
result from the gluon condensate term. Hence, the parameters Sq N ^ and the 
combinations f%(£ N ) straightforwardly extracted from the SR (§) are enhanced 
compared to sg=° « 0.7 GeV 2 and 3/ 2 /(A + 1)(A + 3), resp., by the factors 2 
(for A = 2) and 3 (for N = 4). These are just the results given in Ref. Q. To 
clarify the assumptions implied by such a procedure, we rewrite the CZ sum 
rule using the standard numerical values for the condensates: 
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(A + l)(A + 3) 



0.6 \ 2 „ _ / 4A\ /0.6 X ! 



O.U^ +0.22 1 + - 



(10) 



Taking first A = 0, we see that for M 2 = 0.6 GeV 2 the condensate corrections 
are by factor 3 smaller than the perturbative term while the exponential e~ s l M 
suppresses the A\ contribution by factor 14 compared to the pion one. Hence, 
the sum rule looks very reliable since power corrections are small in the region 
where the s-integral is dominated by the pion. Taking the "first resonance plus 
effective continuum" model for the spectrum and fitting the sum rule in the 
M 2 > 0.6 GeV 2 region gives s ~ 0.75 GeV 2 for the effective threshold, i.e. 
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at the threshold the exponential e~ So//M2 provides 1/3 suppression factor for 
M 2 = 0.6 GeV 2 , which ensures that the result for f 2 is not very sensitive to 
the model chosen for the higher states. 

Now, taking N — 2, we observe that for M 2 — 0.6 GeV 2 the conden- 
sate corrections are by factor 2.4 larger than the perturbative term: the 1/M 2 
expansion is apparently useless at such a value of M 2 . To bring the size of con- 
densate corrections to less than 1/3 of the perturbative term, one should take 
M 2 > 1.2 GeV 2 . However, for such large M 2 values the exponential e~ s l M 
gives practically no suppression at the "old" effective threshold, and results for 
(£ 2 ) would strongly depend on the model for higher states. In particular, the 
"first resonance plus effective continuum" ansatz gives w 1.5 GeV 2 and 
(£ 2 ) « 0.4 which means that with respect to (£ 2 ) the pion is dual to much 
wider interval < s < 1.5 GeV 2 . For N = 4 the duality interval obtained in 
this way is even wider: Sq 4 ' « 2.2 GeV 2 , i.e., the effective continuum threshold 
is assumed to be well above the A\ location. 

Of course, one cannot exclude a priori that a different correlator has a 
different shape of spectral density. Ideally, having the full expression for the 
right-hand side of the sum rule one could find out pn(s) exactly. Having just 
few terms of the 1 /M 2 -expansion, we can only construct an approximation for 
the spectrum, the precision of which depends on the relative magnitude of the 
neglected higher terms. The CZ-proccdure is equivalent to assumption that 
two condensate terms included in their sum rule dominate the expansion not 
only for N = but also for N = 2. In fact, it is impossible to check by a 
direct calculation whether this assumption is true or not, because the number 
of possible condensates explodes when their dimension increases, and there is 
no reliable way to determine their values. Still, it is rather easy to establish 
that coefficients accompanying the condensates (qD 2 q qq) with two covariant 
derivatives D behave like N 3 for large N, i.e. have even larger TV-dependent 
enhancement compared to the perturbative term. In general, the coefficients 
for (q(D 2 ) n qqq) condensates behave like N n+1 for large N. Hence, one would 
rather expect that there are large higher-condensate corrections to the (£ 2 ) 
sum rule. Only some miraculous cancellation can make them small. No reason 
for such a cancellation was given. 

To summarize: if one takes the CZ sum rule at face value, i.e., assumes that 
there are no essential corrections to it, the fitting procedure would produce the 
large CZ value for (£ 2 ). However, since the perturbative term decreases with 
iV while the condensate terms rapidly increase with TV, the CZ sum rules for 
N < 2 is an obvious case when one must expect essential corrections. 
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4 Nonlocal condensates 

It is also instructive to write the SR for the pion DA (pir(x) itself 0: 

+^ 7ra ^-4 g)2 {H[^) + 5(1 - x)] + 2[«5'(x) + - *)]}. (11) 

The OCV) and O(N) terms in Eq. (^|) correspond to the S(x) and (5'(x)-terms 
in Eq. (|l l[) indicating that the vacuum fields are carrying zero fraction of the 
pion momentum. The operator product expansion (underlying eqs.(^|),([l|)) is, 
in fact, a power series expansion over small momenta k of vacuum quarks and 
gluons. Retaining only the (qq) and (GG)-terms (like in eqs.(|g|),(|lT|)) is just 
equivalent to the assumption that k is not simply small but exactly equals zero. 
However, it is much more reasonable to expect that the vacuum quanta have a 
smooth distribution with a finite width /.t. In configuration space, this means 
that vacuum fluctuations have a finite correlation length of the order of 1/fi, so 
that the two-point condensates like (q(0)q(z)) die away for \z\ large compared 
to l//i. In the OPE, (q(Q)q(z)) is expanded in powers of z and the first term 
(g(O)g(O)) produces eventually the <5(x)-term, while higher (q(0)(D 2 ) n q(0)) 
terms give S n (x) contributions resulting in N n factors in the (£ N ) sum rule. In 
other words, arranging the 1/M 2 expansion through the OPE in terms of local 
operators, one automatically obtains (£ N ) in the form of Taylor expansion in N. 
Even if the condensate contribution to the (£ N ) sum rule is a rapidly decreasing 
function of N (which must be the case for any smooth function of £), the OPE 
gives it as a Taylor series in N n whose terms rapidly increase with N. In such a 
situation, it is obviously risky to take just the first term of the expansion, e.g., 
the quark condensate (11 + 47V) factor may well be just the first term of some- 
thing like (11 + 4iV) exp[— N X 2 /M 2 } with much smaller value for N = 2 than 
one would expect from (11 + 47V). How much smaller, depends on the value of 
the scale A 2 . The size of the correlation length of vacuum fluctuations can be 
estimated using the standard valucEj A 2 = (q~D 2 q) / {qq) = 0.4 ± 0.1 GeV 2 for 
the average virtuality of the vacuum quarks. One can see that it is not small 
compared to the relevant hadronic scale Sq =0 w 47r 2 / 2 = 0.7 GeV 2 , and the 
constant-field approximation for the vacuum fields is not safe. Using the expo- 
nential model (q(0)q(z)) — (qq) exp[z 2 A 2 /2] for the nonlocal condensate gives 
a Q£I) sum rule producing the wave functions very close to the asymptotic 
oneO. This study suggests that the humpy form of the CZ wave function 
is a mere consequence of the approximation that vacuum quarks have zero 
momentum. 
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5 QCD sum rule for the 7*7 



form factor 



Another evidence that the pion DA, is close to its asymptotic shape is given by 
a direct QCD sum rule analysis o£3 of the 7*7 — > n° transition form factor. 
In this case, one should consider the three-point correlation function 

^(<Zi,9 2 ) = 2™ J (0\T {j 5 a (Y)J^(X) J„(0)} \Q)e-^ x e lpY d A Xd A Y , (12) 

where J M is the electromagnetic current. The operator product expansion is 
simpler when both photon virtualities q 2 are large: q 2 ,Q 2 > 1 GeV 2 . QCD 
sum rule in this kinematics is given by 



7r/7rF 7 « 7 *7r°(g 2 ; Q 2 ) 



2 I dse 
la 



s/M 2 f 1 xxjxQ 2 +xq 2 )< 



[sxx + xQ 2 + xq 2 ] 3 



dx 



7T 2 ,a. 



GG) 



2M 2 Q 2 2M 2 q 2 Q 2 q< 



64 
243 



w 3 a s (qq} 2 I - 



1 


\Q 2 


9 


9 q 2 ' 


9 9 


M 4 


U 4 


2q 2 


2Q 2 ' Q\ 





(13) 



Keeping only the leading 0(1/ Q 2 and l/g 2 )-terms one can rewrite it as 



F^ {q 2 ,Q 2 ) 

1 



dx 



4tt 



3/ w y (xQ 2 + xq 2 ) 



3M 2 
2tt 2 1 



-8 /M' 



)xx 



2AM 2 x 7T 



— GG)[5{x) +5{x)} 



81M 4 



7ra s (gg) 2 ll[<y(ar) + S(x)} + 2[5'(a;) + S'(x)} 



(14) 



Note, that the expression in curly brackets coincides with the QCD sum rule 
( |lT| ) for the pion DA /^(^(x). Hence, the QCD sum rule approach exactly 
reproduces the PQCD result (|]J). One may be tempted to get a QCD sum rule 
for the integral I by taking q 2 = in Eq. ( |13| ) ■ The attempt is ruined by power 
singularities 1/q 2 , l/g 4 in the condensate terms. Moreover, the perturbative 
term in the small-g 2 region has logarithms log q 2 which are a typical example of 
mass singularities (see, e.g.,E3). All these infrared sensitive terms are produced 
in a regime when the hard momentum flow bypasses the soft photon vertex, 
i.e., the EM current J^X) of the low-virtuality photon is far away from the 
two other currents J(0), j (Y). It is also important to observe that power 
singularities 1/q 2 , 1/q 4 are generated precisely by the same S(x) and S'(x) 
terms in Eq.(|l4|) which generate the two- hump form for in the CZ- 

approachQ. As shown in Reflii, the humps disappear if one treats the S(x) 
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and S'(x) terms as the first terms of a formal expansion $(a;) ~ ^2 &n5 n {x) of 
smooth functions $>(x). Similarly, the 1/q 2 singularity can be understood as 
the first term of the large-g 2 expansion of a term like l/(q 2 + m 2 ) in powers 
of l/q 2 . However, constructing <j>(x) from two first terms of such expansion 
is a strongly model-dependent procedure. On the other hand, the small-q 2 
behavior of the three-point function is rather severely constrained by known 
structure of the physical spectrum in the EM-current channel. The procedure 
developed in RefsOuJ allows to subtract all the small-g 2 singularities from the 
coefficient functions of the original OPE for the 3-point correlation function 
Eq. (|l2|) . They are absorbed in this approach by universal bilocal correlators, 
which can be also interpreted as moments of the DAs for (almost) real photon 



e^ x (0\T{JJX)O^(0)}\0)d 4 X, 



where On ^ (0) acepOperators of leading and next-to-leading twist with n covari- 
ant derivatives oE_3. The bilocal contribution to the 3-point function Eq.(|l2|) 
can be written in a "parton" form as a convolution of the photon DAs and 
some coefficient functions. The latter originate from a light cone OPE for the 
product T{ J(0)j 5 (Y)}. The amplitude T is now a sum of its purely short- 
distance (SD) (regular for q 2 = 0) and bilocal (B) parts. Getting the q 2 — > 

limit of Hn\qi) requires a nonperturbative input obtained from an auxiliary 
QCD sum rule. After all the modifications outlined above are made, one can 
write the QCD sum rule for the 77* — > ir° form factor in the q 2 = limit: 



1 - 2 



Q 2 - 2s 



4 -6 S Q 2 + 3s 2 ( s 2 



9 ' 7T 



{a + Q 2 ) 4 
GG) 



1 



2 3m 2 
1 



2Q 2 M 2 Q 4 
1 



64 c, ._ , 2 ,. 

H 7r ct.(qq) lim . , 

27 wy/ x^a{2Q 2 M 4 Q 4 m 2 p 






4 

0® 



s 2 + 3sQ 2 + 4Q 4 



(s + Q 2 ) 3 



( Sp 2 


m 2 J 


i/M 2 Q 


) 2 ds 




-Q 2 ) 2 


-s/M 2 


ds 




+ Q 2 ) 
- 2 


x2 )ds 





log — - 3 



^ _ s/M 2 ^ 2 + 3sQ 2 + 6Q 4 

2\3 



(s + Q 2 ) 



(15) 



+ A 2 
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Here the bilocal contributions are modeled by asymptotic form for the p-meson 
DAs. They are approximately dual to the_cprresponding perturbative contri- 
bution with the p- meson duality interval E3 s p = 1.5 GeV 2 . The results of 
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Figure 1: 

fitting prnr-prliirp for ( [151) favor the value so ~ 0.7 GeV 2 as the effective thresh- 
oldoEj. For this reason, the results of our calculations are well approximated 
by the local duality prescriptionO 



ttUF^o(Q 2 ) 



quark 



(s,Q 2 )ds = 



1 



l + Q 2 /s 



(16) 



which coincides for sq — 47r 2 / 2 with the BL interpolation formula. In Fig.[l], we 
present our curve (solid line) for Q 2 F 11 * 7r o(Q 2 )/A-Kf. K calculated from Eq.(p"5|) 
for so = 0.7 GeV 2 and M 2 — OJ&GeV 2 . One can observe very good agree- 
ment with the new CLEO data El It is also rather close to the BL interpo- 
lation/local duality formula (long-dashed line) and the p-pole approximation 
(short-dashed line) nUF VMD (Q 2 ) = 1/(1 + Q 2 /m 2 p ). It should be noted that 
the <5 2 -dependence of the p-pole type emerges due to the fact that the pion 
duality interval so ~ 0.67 GeV 2 is numerically close to m 2 w 0.6 GeV 2 . In 
the region Q 2 > Q 2 ~ 3 GeV 2 , our curve for Q 2 F 1 ^, w o(Q 2 ) is practically con- 
stant, supporting the PQCD expectation (Q). The absolute magnitude of our 
prediction gives / w 2.4 for the /-integral with an accuracy of about 20%. Com- 
paring the value I = 2.4 with J as = 3 and I cz = 5, we conclude that our result 
favours a pion DA which is narrower than the asymptotic form. Parametriz- 
ing the width of (p^(x) by a simple model ip n (x) ~ [a;(l — x)] n , we obtain that 



/ = 2.4 corresponds to n — 2.5. The second moment (£ 2 ) = ((a; 



'} for such 
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a function is 0.125 (recall that (£ 2 } as = 0.2 while (£ 2 ) cz = 0-43) which agrees 
with the lattice calculation E3. 

Thus, the old claim EJ that the CZ sum rules □ for the moments of DAs 
are unreliable is now sunnoiled both by a direct QCD sum rule calculation 
of the 7*77T° form factorOcJ producing the result corresponding—to a narrow 
pion DA, and by experimental measurement of this form factor 13 which also 
favors a pion DA close to the asymptotic form. Since the humpy form of the 
CZ models for the nucleon DA'sEl has the same origin as in the pion case, there 
is no doubt that the nucleon DA's are also close to the asymptotic ones. This 
means that PQCD contributions to nucleon elastic and transition form factors 
are tiny at available and reachable energies. 
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